A Flexible Method for Identifying Spatial Clusters of Breast Cancer Using Individual-Level Data.

Prior research has shown that cancer risk varies by geography, but scan statistics methods for identifying cancer clusters in case-control studies have been limited in their ability to identify multiple clusters and adjust for participant-level risk factors. We develop a method to identify geographic patterns of breast cancer odds using the Wisconsin Women's Health Study, a series of 5 population-based case-control studies of female Wisconsin residents aged 20-79 enrolled in 1988-2004 (cases=16,076, controls=16,795). We create sets of potential clusters by overlaying a 1 km grid over each county-neighborhood and enumerating a series of overlapping circles. Using a two-step approach, we fit a penalized binomial regression model to the number of cases and trials in each grid cell, penalizing all potential clusters by the least absolute shrinkage and selection operator (Lasso). We use BIC to select the number of clusters, which are included in a participant-level logistic regression model. We identify 15 geographic clusters, resulting in 23 areas of unique geographic odds ratios. After adjustment for known risk factors, confidence intervals narrowed but breast cancer odds ratios did not meaningfully change; one additional hotspot was identified. By considering multiple overlapping spatial clusters simultaneously, we discern gradients of spatial odds across Wisconsin.

Regularized spatial and spatio-temporal cluster detection.

Spatial and spatio-temporal cluster detection are important tools in public health and many other areas of application. Cluster detection can be approached as a multiple testing problem, typically using a space and time scan statistic. We recast the spatial and spatio-temporal cluster detection problem in a high-dimensional data analytical framework with Poisson or quasi-Poisson regression with the Lasso penalty. We adopt a fast and computationally-efficient method using a novel sparse matrix representation of the effects of potential clusters. The number of clusters and tuning parameters are selected based on (quasi-)information criteria. We evaluate the performance of our proposed method including the false positive detection rate and power using a simulation study. Application of the method is illustrated using breast cancer incidence data from three prefectures in Japan.

The number of tree species on Earth.

One of the most fundamental questions in ecology is how many species inhabit the Earth. However, due to massive logistical and financial challenges and taxonomic difficulties connected to the species concept definition, the global numbers of species, including those of important and well-studied life forms such as trees, still remain largely unknown. Here, based on global ground-sourced data, we estimate the total tree species richness at global, continental, and biome levels. Our results indicate that there are ∼73,000 tree species globally, among which ∼9,000 tree species are yet to be discovered. Roughly 40% of undiscovered tree species are in South America. Moreover, almost one-third of all tree species to be discovered may be rare, with very low populations and limited spatial distribution (likely in remote tropical lowlands and mountains). These findings highlight the vulnerability of global forest biodiversity to anthropogenic changes in land use and climate, which disproportionately threaten rare species and thus, global tree richness.

Clustered spatio-temporal varying coefficient regression model.

In regression analysis for spatio-temporal data, identifying clusters of spatial units over time in a regression coefficient could provide insight into the unique relationship between a response and covariates in certain subdomains of space and time windows relative to the background in other parts of the spatial domain and the time period of interest. In this article, we propose a varying coefficient regression method for spatial data repeatedly sampled over time, with heterogeneity in regression coefficients across both space and over time. In particular, we extend a varying coefficient regression model for spatial-only data to spatio-temporal data with flexible temporal patterns. We consider the detection of a potential cylindrical cluster of regression coefficients based on testing whether the regression coefficient is the same or not over the entire spatial domain for each time point. For multiple clusters, we develop a sequential identification approach. We assess the power and identification of known clusters via a simulation study. Our proposed methodology is illustrated by the analysis of a cancer mortality dataset in the Southeast of the U.S.

Subjective sleep measurement: comparing sleep diary to questionnaire.

PURPOSE: The sleep diary is the gold standard of self-reported sleep duration, but its comparability to sleep questionnaires is uncertain. The purpose of this study was to compare self-reported sleep duration between a sleep diary and a sleep questionnaire and to test whether sleep-related disorders were associated with diary-questionnaire differences in sleep duration. PARTICIPANTS AND METHODS: We compared self-reported sleep duration from 5,432 questionnaire-sleep diary pairs in a longitudinal cohort of 1,516 adults. Participants reported sleep information in seven-day sleep diaries and in questionnaires. Research staff abstracted average sleep durations for three time periods (overall; weekday; weekend) from diaries and questionnaires. For each time period, we evaluated diary-questionnaire differences in sleep duration with Welch's two-sample t-tests. Using linear mixed effects regression, we regressed overall diary-questionnaire sleep duration difference on several participant characteristics: reporting any insomnia symptoms, having sleep apnea, sex, body mass index, smoking status, Short Form-12 Physical Health Composite Score, and Short Form-12 Mental Health Composite Score. RESULTS: The average diary-reported overall sleep duration (7.76 hrs) was longer than that of the questionnaire (7.07 hrs) by approximately 41 mins (0.69 hrs, 95% confidence interval: 0.62, 0.76 hrs). Results were consistent across weekday- and weekend-specific differences. Demographic-adjusted linear mixed effects models tested whether insomnia symptoms or sleep apnea were associated with diary-questionnaire differences in sleep duration. Insomnia symptoms were associated with a 17 min longer duration on the diary relative to the questionnaire (ß=0.28 hrs, 95% confidence interval: 0.22, 0.33 hrs), but sleep apnea was not significantly associated with diary-questionnaire difference. Female sex was associated with greater diary-questionnaire duration differences, whereas better self-reported health was associated with lesser differences. CONCLUSION: Diaries and questionnaires are somewhat disparate methods of assessing subjective sleep duration, although diaries report longer duration relative to questionnaires, and insomnia symptoms may contribute to greater perceived differences.

Spatial Regression Analysis of Poverty in R.

Poverty has been studied across many social science disciplines, resulting in a large body of literature. Scholars of poverty research have long recognized that the poor are not uniformly distributed across space. Understanding the spatial aspect of poverty is important because it helps us understand place-based structural inequalities. There are many spatial regression models, but there is a learning curve to learn and apply them to poverty research. This manuscript aims to introduce the concepts of spatial regression modeling and walk the reader through the steps of conducting poverty research using R: standard exploratory data analysis, standard linear regression, neighborhood structure and spatial weight matrix, exploratory spatial data analysis, and spatial linear regression. We also discuss the spatial heterogeneity and spatial panel aspects of poverty. We provide code for data analysis in the R environment and readers can modify it for their own data analyses. We also present results in their raw format to help readers become familiar with the R environment.

The ruminal bacterial community in lactating dairy cows has limited variation on a day-to-day basis.

Dairy cows rely on a complex ruminal microbiota to digest their host-indigestible feed. Our ability to characterize this microbiota has advanced significantly due to developments in next-generation sequencing. However, efforts to sample the rumen, which typically involves removing digesta directly from the rumen via a cannula, intubation, or rumenocentesis, is costly and labor intensive. As a result, the majority of studies characterizing the rumen microbiota are conducted on samples collected at a single time point. Currently, it is unknown whether there is significant day-to-day variation in the rumen microbiota, a factor that could strongly influence conclusion drawn from studies that sample at a single time point. To address this, we examined day-to-day changes in the ruminal microbiota of lactating dairy cows using next-generation sequencing to determine if single-day sampling is representative of sampling across 3 consecutive days. We sequenced single-day solid and liquid fractions of ruminal digesta collected over 3 consecutive days from 12 cannulated dairy cows during the early, middle, and late stages of a single lactation cycle using the V4 region of the bacterial 16S rRNA gene. We then generated 97% similarity operational taxonomic units (OTUs) from these sequences and showed that any of the individual samples from a given 3-day sampling period is equivalent to the mean OTUs determined from the combined 3-d data set. This finding was consistent for both solid and liquid fractions of the rumen, and we thus conclude that there is limited day-to-day variability in the rumen microbiota.

Control of phosphorothioate stereochemistry substantially increases the efficacy of antisense oligonucleotides.

Whereas stereochemical purity in drugs has become the standard for small molecules, stereoisomeric mixtures containing as many as a half million components persist in antisense oligonucleotide (ASO) therapeutics because it has been feasible neither to separate the individual stereoisomers, nor to synthesize stereochemically pure ASOs. Here we report the development of a scalable synthetic process that yields therapeutic ASOs having high stereochemical and chemical purity. Using this method, we synthesized rationally designed stereopure components of mipomersen, a drug comprising 524,288 stereoisomers. We demonstrate that phosphorothioate (PS) stereochemistry substantially affects the pharmacologic properties of ASOs. We report that Sp-configured PS linkages are stabilized relative to Rp, providing stereochemical protection from pharmacologic inactivation of the drug. Further, we elucidated a triplet stereochemical code in the stereopure ASOs, 3'-SpSpRp, that promotes target RNA cleavage by RNase H1 in vitro and provides a more durable response in mice than stereorandom ASOs.

Effect of an RNA interference drug on the synthesis of proprotein convertase subtilisin/kexin type 9 (PCSK9) and the concentration of serum LDL cholesterol in healthy volunteers: a randomised, single-blind, placebo-controlled, phase 1 trial.

BACKGROUND: Proprotein convertase subtilisin/kexin type 9 (PCSK9) binds to LDL receptors, leading to their degradation. Genetics studies have shown that loss-of-function mutations in PCSK9 result in reduced plasma LDL cholesterol and decreased risk of coronary heart disease. We aimed to investigate the safety and efficacy of ALN-PCS, a small interfering RNA that inhibits PCSK9 synthesis, in healthy volunteers with raised cholesterol who were not on lipid-lowering treatment. METHODS: We did a randomised, single-blind, placebo-controlled, phase 1 dose-escalation study in healthy adult volunteers with serum LDL cholesterol of 3·00 mmol/L or higher. Participants were randomly assigned in a 3:1 ratio by computer algorithm to receive one dose of intravenous ALN-PCS (with doses ranging from 0·015 to 0·400 mg/kg) or placebo. The primary endpoint was safety and tolerability of ALN-PCS. Secondary endpoints were the pharmacokinetic characteristics of ALN-PCS and its pharmacodynamic effects on PCSK9 and LDL cholesterol. Study participants were masked to treatment assignment. Analysis was per protocol and we used ANCOVA to analyse pharmacodynamic endpoint data. This trial is registered with ClinicalTrials.gov, number NCT01437059. FINDINGS: Of 32 participants, 24 were randomly allocated to receive a single dose of ALN-PCS (0·015 mg/kg [n=3], 0·045 mg/kg [n=3], 0·090 mg/kg [n=3], 0·150 mg/kg [n=3], 0·250 mg/kg [n=6], or 0·400 mg/kg [n=6]) and eight to placebo. The proportions of patients affected by treatment-emergent adverse events were similar in the ALN-PCS and placebo groups (19 [79%] vs seven [88%]). ALN-PCS was rapidly distributed, with peak concentration and area under the curve (0 to last measurement) increasing in a roughly dose-proportional way across the dose range tested. In the group given 0·400 mg/kg of ALN-PCS, treatment resulted in a mean 70% reduction in circulating PCSK9 plasma protein (p<0·0001) and a mean 40% reduction in LDL cholesterol from baseline relative to placebo (p<0·0001). INTERPRETATION: Our results suggest that inhibition of PCSK9 synthesis by RNA interference (RNAi) provides a potentially safe mechanism to reduce LDL cholesterol concentration in healthy individuals with raised cholesterol. These results support the further assessment of ALN-PCS in patients with hypercholesterolaemia, including those being treated with statins. This study is the first to show an RNAi drug being used to affect a clinically validated endpoint (ie, LDL cholesterol) in human beings. FUNDING: Alnylam Pharmaceuticals.

Silencing of lipid metabolism genes through IRE1α-mediated mRNA decay lowers plasma lipids in mice.

XBP1 is a key regulator of the unfolded protein response (UPR), which is involved in a wide range of physiological and pathological processes. XBP1 ablation in liver causes profound hypolipidemia in mice, highlighting its critical role in lipid metabolism. XBP1 deficiency triggers feedback activation of its upstream enzyme IRE1α, instigating regulated IRE1-dependent decay (RIDD) of cytosolic mRNAs. Here, we identify RIDD as a crucial control mechanism of lipid homeostasis. Suppression of RIDD by RNA interference or genetic ablation of IRE1α reversed hypolipidemia in XBP1-deficient mice. Comprehensive microarray analysis of XBP1 and/or IRE1α-deficient liver identified genes involved in lipogenesis and lipoprotein metabolism as RIDD substrates, which might contribute to the suppression of plasma lipid levels by activated IRE1α. Ablation of XBP1 ameliorated hepatosteatosis, liver damage, and hypercholesterolemia in dyslipidemic animal models, suggesting that direct targeting of either IRE1α or XBP1 might be a feasible strategy to treat dyslipidemias.

Activation of ER stress and mTORC1 suppresses hepatic sortilin-1 levels in obese mice.

Recent GWAS have identified SNPs at a human chromosom1 locus associated with coronary artery disease risk and LDL cholesterol levels. The SNPs are also associated with altered expression of hepatic sortilin-1 (SORT1), which encodes a protein thought to be involved in apoB trafficking and degradation. Here, we investigated the regulation of Sort1 expression in mouse models of obesity. Sort1 expression was markedly repressed in both genetic (ob/ob) and high-fat diet models of obesity; restoration of hepatic sortilin-1 levels resulted in reduced triglyceride and apoB secretion. Mouse models of obesity also exhibit increased hepatic activity of mammalian target of rapamycin complex 1 (mTORC1) and ER stress, and we found that administration of the mTOR inhibitor rapamycin to ob/ob mice reduced ER stress and increased hepatic sortilin-1 levels. Conversely, genetically increased hepatic mTORC1 activity was associated with repressed Sort1 and increased apoB secretion. Treating WT mice with the ER stressor tunicamycin led to marked repression of hepatic sortilin-1 expression, while administration of the chemical chaperone PBA to ob/ob mice led to amelioration of ER stress, increased sortilin-1 expression, and reduced apoB and triglyceride secretion. Moreover, the ER stress target Atf3 acted at the SORT1 promoter region as a transcriptional repressor, whereas knockdown of Atf3 mRNA in ob/ob mice led to increased hepatic sortilin-1 levels and decreased apoB and triglyceride secretion. Thus, in mouse models of obesity, induction of mTORC1 and ER stress led to repression of hepatic Sort1 and increased VLDL secretion via Atf3. This pathway may contribute to dyslipidemia in metabolic disease.

The Scap/SREBP pathway is essential for developing diabetic fatty liver and carbohydrate-induced hypertriglyceridemia in animals.

Insulin resistance leads to hypertriglyceridemia and hepatic steatosis and is associated with increased SREBP-1c, a transcription factor that activates fatty acid synthesis. Here, we show that steatosis in insulin-resistant ob/ob mice was abolished by deletion of Scap, an escort protein necessary for generating nuclear isoforms of all three SREBPs. Scap deletion reduced lipid synthesis and prevented fatty livers despite persistent obesity, hyperinsulinemia, and hyperglycemia. Scap deficiency also prevented steatosis in mice fed high-fat diets. Steatosis was also prevented when siRNAs were used to silence Scap in livers of sucrose-fed hamsters, a model of diet-induced steatosis and hypertriglyceridemia. This silencing reduced all three nuclear SREBPs, decreasing lipid biosynthesis and abolishing sucrose-induced hypertriglyceridemia. These results demonstrate that SREBP activation is essential for development of diabetic hepatic steatosis and carbohydrate-induced hypertriglyceridemia, but not insulin resistance. Inhibition of SREBP activation has therapeutic potential for treatment of hypertriglyceridemia and fatty liver disease.

IRE1α activation protects mice against acetaminophen-induced hepatotoxicity.

The mammalian stress sensor IRE1α plays a central role in the unfolded protein, or endoplasmic reticulum (ER), stress response by activating its downstream transcription factor XBP1 via an unconventional splicing mechanism. IRE1α can also induce the degradation of a subset of mRNAs in a process termed regulated IRE1-dependent decay (RIDD). Although diverse mRNA species can be degraded by IRE1α in vitro, the pathophysiological functions of RIDD are only beginning to be explored. Acetaminophen (APAP) overdose is the most frequent cause of acute liver failure in young adults in the United States and is primarily caused by CYP1A2-, CYP2E1-, and CYP3A4-driven conversion of APAP into hepatotoxic metabolites. We demonstrate here that genetic ablation of XBP1 results in constitutive IRE1α activation in the liver, leading to RIDD of Cyp1a2 and Cyp2e1 mRNAs, reduced JNK activation, and protection of mice from APAP-induced hepatotoxicity. A pharmacological ER stress inducer that activated IRE1α suppressed the expression of Cyp1a2 and Cyp2e1 in WT, but not IRE1α-deficient mouse liver, indicating the essential role of IRE1α in the down-regulation of these mRNAs upon ER stress. Our study reveals an unexpected function of RIDD in drug metabolism.

From noncoding variant to phenotype via SORT1 at the 1p13 cholesterol locus.

Recent genome-wide association studies (GWASs) have identified a locus on chromosome 1p13 strongly associated with both plasma low-density lipoprotein cholesterol (LDL-C) and myocardial infarction (MI) in humans. Here we show through a series of studies in human cohorts and human-derived hepatocytes that a common noncoding polymorphism at the 1p13 locus, rs12740374, creates a C/EBP (CCAAT/enhancer binding protein) transcription factor binding site and alters the hepatic expression of the SORT1 gene. With small interfering RNA (siRNA) knockdown and viral overexpression in mouse liver, we demonstrate that Sort1 alters plasma LDL-C and very low-density lipoprotein (VLDL) particle levels by modulating hepatic VLDL secretion. Thus, we provide functional evidence for a novel regulatory pathway for lipoprotein metabolism and suggest that modulation of this pathway may alter risk for MI in humans. We also demonstrate that common noncoding DNA variants identified by GWASs can directly contribute to clinical phenotypes.

Targeted delivery of RNAi therapeutics with endogenous and exogenous ligand-based mechanisms.

Lipid nanoparticles (LNPs) have proven to be highly efficient carriers of short-interfering RNAs (siRNAs) to hepatocytes in vivo; however, the precise mechanism by which this efficient delivery occurs has yet to be elucidated. We found that apolipoprotein E (apoE), which plays a major role in the clearance and hepatocellular uptake of physiological lipoproteins, also acts as an endogenous targeting ligand for ionizable LNPs (iLNPs), but not cationic LNPs (cLNPs). The role of apoE was investigated using both in vitro studies employing recombinant apoE and in vivo studies in wild-type and apoE(-/-) mice. Receptor dependence was explored in vitro and in vivo using low-density lipoprotein receptor (LDLR(-/-))-deficient mice. As an alternative to endogenous apoE-based targeting, we developed a targeting approach using an exogenous ligand containing a multivalent N-acetylgalactosamine (GalNAc)-cluster, which binds with high affinity to the asialoglycoprotein receptor (ASGPR) expressed on hepatocytes. Both apoE-based endogenous and GalNAc-based exogenous targeting appear to be highly effective strategies for the delivery of iLNPs to liver.

Lipid-like materials for low-dose, in vivo gene silencing.

Significant effort has been applied to discover and develop vehicles which can guide small interfering RNAs (siRNA) through the many barriers guarding the interior of target cells. While studies have demonstrated the potential of gene silencing in vivo, improvements in delivery efficacy are required to fulfill the broadest potential of RNA interference therapeutics. Through the combinatorial synthesis and screening of a different class of materials, a formulation has been identified that enables siRNA-directed liver gene silencing in mice at doses below 0.01 mg/kg. This formulation was also shown to specifically inhibit expression of five hepatic genes simultaneously, after a single injection. The potential of this formulation was further validated in nonhuman primates, where high levels of knockdown of the clinically relevant gene transthyretin was observed at doses as low as 0.03 mg/kg. To our knowledge, this formulation facilitates gene silencing at orders-of-magnitude lower doses than required by any previously described siRNA liver delivery system.

Therapeutic RNAi targeting PCSK9 acutely lowers plasma cholesterol in rodents and LDL cholesterol in nonhuman primates.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) regulates low density lipoprotein receptor (LDLR) protein levels and function. Loss of PCSK9 increases LDLR levels in liver and reduces plasma LDL cholesterol (LDLc), whereas excess PCSK9 activity decreases liver LDLR levels and increases plasma LDLc. Here, we have developed active, cross-species, small interfering RNAs (siRNAs) capable of targeting murine, rat, nonhuman primate (NHP), and human PCSK9. For in vivo studies, PCSK9 and control siRNAs were formulated in a lipidoid nanoparticle (LNP). Liver-specific siRNA silencing of PCSK9 in mice and rats reduced PCSK9 mRNA levels by 50-70%. The reduction in PCSK9 transcript was associated with up to a 60% reduction in plasma cholesterol concentrations. These effects were shown to be mediated by an RNAi mechanism, using 5'-RACE. In transgenic mice expressing human PCSK9, siRNAs silenced the human PCSK9 transcript by >70% and significantly reduced PCSK9 plasma protein levels. In NHP, a single dose of siRNA targeting PCSK9 resulted in a rapid, durable, and reversible lowering of plasma PCSK9, apolipoprotein B, and LDLc, without measurable effects on either HDL cholesterol (HDLc) or triglycerides (TGs). The effects of PCSK9 silencing lasted for 3 weeks after a single bolus i.v. administration. These results validate PCSK9 targeting with RNAi therapeutics as an approach to specifically lower LDLc, paving the way for the development of PCSK9-lowering agents as a future strategy for treatment of hypercholesterolemia.

Gene silencing activity of siRNAs with a ribo-difluorotoluyl nucleotide.

Recently, chemically synthesized short interfering RNA (siRNA) duplexes have been used with success for gene silencing. Chemical modification is desired for therapeutic applications to improve biostability and pharmacokinetic properties; chemical modification may also provide insight into the mechanism of silencing. siRNA duplexes containing the 2,4-difluorotoluyl ribonucleoside (rF) were synthesized to evaluate the effect of noncanonical nucleoside mimetics on RNA interference. 5'-Modification of the guide strand with rF did not alter silencing relative to unmodified control. Internal uridine to rF substitutions were well-tolerated. Thermal melting analysis showed that the base pair between rF and adenosine (A) was destabilizing relative to a uridine-adenosine pair, although it was slightly less destabilizing than other mismatches. The crystal structure of a duplex containing rFoA pairs showed local structural variations relative to a canonical RNA helix. As the fluorine atoms cannot act as hydrogen bond acceptors and are more hydrophobic than uridine, there was an absence of a well-ordered water structure around the rF residues in both grooves. siRNAs with the rF modification effectively silenced gene expression and offered improved nuclease resistance in serum; therefore, evaluation of this modification in therapeutic siRNAs is warranted.

Small-molecule modulators of Hedgehog signaling: identification and characterization of Smoothened agonists and antagonists.

BACKGROUND: The Hedgehog (Hh) signaling pathway is vital to animal development as it mediates the differentiation of multiple cell types during embryogenesis. In adults, Hh signaling can be activated to facilitate tissue maintenance and repair. Moreover, stimulation of the Hh pathway has shown therapeutic efficacy in models of neuropathy. The underlying mechanisms of Hh signal transduction remain obscure, however: little is known about the communication between the pathway suppressor Patched (Ptc), a multipass transmembrane protein that directly binds Hh, and the pathway activator Smoothened (Smo), a protein that is related to G-protein-coupled receptors and is capable of constitutive activation in the absence of Ptc. RESULTS: We have identified and characterized a synthetic non-peptidyl small molecule, Hh-Ag, that acts as an agonist of the Hh pathway. This Hh agonist promotes cell-type-specific proliferation and concentration-dependent differentiation in vitro, while in utero it rescues aspects of the Hh-signaling defect in Sonic hedgehog-null, but not Smo-null, mouse embryos. Biochemical studies with Hh-Ag, the Hh-signaling antagonist cyclopamine, and a novel Hh-signaling inhibitor Cur61414, reveal that the action of all these compounds is independent of Hh-protein ligand and of the Hh receptor Ptc, as each binds directly to Smo. CONCLUSIONS: Smo can have its activity modulated directly by synthetic small molecules. These studies raise the possibility that Hh signaling may be regulated by endogenous small molecules in vivo and provide potent compounds with which to test the therapeutic value of activating the Hh-signaling pathway in the treatment of traumatic and chronic degenerative conditions.